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<3 ABSTRACT 

The H.E.S.S. very-high-energy (VHE, E > 0.1 TeV) gamma-ray telescope system has discovered a new source, HESS J1747-248. The measured 
integral flux is (1.2 ± 0.3) x 10~ 12 cnT 2 s~' above 440 GeV for a power-law photon spectral index of 2.5 ± 0.3 stat ± 0.2 sys . The VHE gamma-ray 
source is located in the close vicinity of the Galactic globular cluster Terzan 5 and extends beyond the H.E.S.S. point spread function (0.07°). The 
probability of a chance coincidence with Terzan 5 and an unrelated VHE source is quite low (~ 10~ 4 ). With the largest population of identified 
millisecond pulsars (msPSRs), a very high core stellar density and the brightest GeV range flux as measured by Fermi-LXT, Terzan 5 stands out 
■ among Galactic globular clusters. The properties of the VHE source are briefly discussed in the context of potential emission mechanisms, notably 
- 1 in relation to msPSRs. Interpretation of the available data accommodates several possible origins for this VHE gamma-ray source, although none 
of them offers a satisfying explanation of its peculiar morphology. 

Key words. Galaxy:globular clusters: individual: Terzan 5 - Radiation mechanisms: non-thermal - (Stars:) pulsars: general - Gamma rays: 
general 
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1. Introduction 

Several types of Galactic VHE y-ray sources, such as pulsar 
wind nebulae (PWNe) and supernova remnants, have been de- 
tected to date, but so far none are in the vicinity of a glob- 
ular cluster (GC). GCs are very old stellar systems with ex- 
ceptionally high densities of stars in their cores, leading to nu- 
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merous stellar collisions (see e.g. lPoolev & Hut 2006 ), and the y 
also contain many millisecond pulsars (msPSRs, lRansoml2008l) . 
likely related to the la rge number of binary stellar members 
dCamilo & Rasioll2005l) . GCs are predicted to emit VHE y-rays 
from inverse Compton (IC) up-scattering of photons from both 
stellar radiation fields and the cosmic microwave background by 
energetic electrons. Electrons could be produced at least at two 
differ ent sites: in the magnetosphere of the msPSRs (IVenter et al.l 
[20091 hereafter VJ09) or at msPSR wind nebulae shocks where 
electrons produced by the msPSRs dBednarek & Sitarek1l2007l 
hereafter BS07) could be re-accelerated. Recently observed 
GeV y -ray emission from several GCs (Fermi-LAT observa- 
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iTam et af[ 201 lb prove that electrons accelerated up to the re 



quired energy range are present in GCs. 

This paper reports the discovery of VHE y-rays from the di- 
rection of the GC Terzan 5. Terzan 5 hosts t he largest popu la- 
tion of msPSRs detected in a GC so far (33, lRansomll2008b . It 
also contains the largest estimated number of msPSRs as de- 
rived from its emission in th e GeV ( 1 80 + '2°. lAbdo et al.ll2010al) 
and ra dio domains (60 - 200. lFruchter & Gossl2000t Kong et all 
120101) . It is located at a distance of 5.9 kpc dFerraro et al.ll2009h 
at RA(J2000) 17 h 48 m 04 s .85 and Dec -24°46'44".6 (Galactic co- 
ordinates: / = 3.84°, b - 1.69°) and exhibits a core radius 
r c = 0'.15 , a half-mass radius r^ = 0'.52, and a tidal radius 
r t = 4'.6 dLanzoni et al. I I20T0I) . The detection of diffuse non- 
thermal X-ray emission centred on the GC bu t extending beyon d 
rh (between 1 and 5.8 rh)has been reported bv lEger et al.l(l2010l) . 
Observations of other GCs in the VHE domain so far have only 
resulted in upper hmits (on M 13, M 15, M 5, and 47 Tucanae in 
lAnderhub et al.l2009LlMcCutcheonl2009l:lAharonian et al.l2009l 
respectively, and references therein for older results). 



2. Observation and analysis 

H.E.S.S. is an array of four Imaging Atmospheric Cherenkov 
Telescopes, located in the Khomas Highland of Namibia (lat- 
itude 23.16°South, altitude 1800 m). Stereoscopic trigger and 
analysis methods allow efficient background (cosmic ray, CR) 
rejection and accurate reconstruction of energy (better than 20%) 
and arrival direction (better than 0. l°per event) for y-ray s in the 
range 0. 1 - 100 TeV For point-like sources, the array has a nom- 
inal detection sensitivity of ~ 1 % of the flux of the Crab Nebula 
(Crab) above 1 TeV with a significance of 5 cr in 25 hours of 
observation at small zenith angles (< 20°). A thorough dis- 
cussion of the H.E.S.S. standard analysis based on Hillas pa- 
ra meters and the perform ance of the instrument can be found 
in lAharonian et al.l d2006l) . Data on Terzan 5 was obtained by 
H.E.S.S. from 2004 to 2010 both as part of its systematic survey 
of the Galactic plane and with dedicated observations. 

For source detection, the data quality cuts result in 90 hours 
of 3- and 4-telescope data with an average zenith angle of 
20.4°and a mean pointing directio n offset from Terzan 5 of 
0.95°. After applying hard cuts (see lAharonian et ai1l2006l) op- 
timized for point-like sources, with a corresponding energy 
threshold of 380 GeV and a point spread function (PSF) 68% 
containment radius Tpsf = 0.07° (4'. 2), a source of VHE y-ray s 
is detected (see Fig.Q] produce d using the tem plate background 
estimation method, described in lRowellll200l . The significance 
reaches 5.3 cr at the position of Terzan 5 with a nearby peak 
significance of 7.5 cr pre-trial (above 5 cr post-trial). All results 
are confirmed by an alysis based on fully independ ent calibration 
and analysis chains dde Naurois & Rollandll2009l) . 
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Fig. 1. Exposure-corrected excess image from the H.E.S.S. data, 
smoothed with a Gaussian function of width 0.1°and overlaid 
with significance contours (4 - 6 cr) in RADec J2000 coordi- 
nates. The circles show the half-mass radius (in black) and the 
larger tidal radius (in cyan) of the GC. The cross indicates the 
best-fit source position of HESS J1747 - 248, assuming a 2D 
Gaussian shape, with 1 cr uncertainty on each axis. The rect- 
angle represents the integration region used for the full-source 
spectral analysis. The upper-right corner circle illustrates the in- 
strumental PSF. 



The source appears to extend beyond the GC tidal radius. A 
2D Gaussian fit 1 results in a best-fit position (^ 2 /ndf = 0.18) 
RA(J2000) 17 h 47 m 4 9 s ± l m .8 jya ,± V3 SVS and Dec -24°48'30" + 
36" stat ± 20" sys (see lAcero et"alll2010al) . offset by 4'.0 + 1'.9 
from the GC centre, therefore, this new VHE y-ray source is 
named HESS J1747 - 248. The size of the source is given by 
the Gaussian widths 9'. 6 + 2'.4 and 1'.8 ± 1'.2 for the major and 
minor axes, respectively, oriented 92° ±6° westwards from north. 

To establish the photon spectrum of the source, a more re- 
strictive data selection is applied, to improve the energy recon- 
struction, resulting in a total of 62 h of live time. Figure [2] illus- 
trates the results of a spect ral analysis with th e reflected back- 
ground estimation method dBerge et alj|2007l) . The test region 
(shown in Fig. [TJ is a rectangle defined from the 2D Gaussian fit 
(half dimensions rpsF + cr x , 2 x rpsF, rotated and centred on the 
best-fit position). Larger integration regions give lower signal 

significance. For a power-law spectral model ^(^) , the flux 

normalization k at 1 TeV is (5.2+1. I)xl0~ 13 cm" 2 s" 1 TeV -1 , 
and the spectral index F = 2.5 + 0.3 stat ± 0.2 sys , correspond- 
ing to an integral photon flux within the integration region of 
(1 .2+0.3)x 10 -12 cnr 2 s -1 , or 1 .5% of the Crab flux, in the 0.44 - 
24 TeV range (^- 2 /ndf = 1.1). There are not enough excess 
events to discuss a more complex spectral model. 

For comparis on with the VHE y-ray flux upper limit on the 
GC 47 Tucanae dAharonian et al.l 120091) . see Sect. [33] where 
a spectral analysis for a point-like source centred on the core 
of Terzan 5 is carried out. It results in a compatible photon 
index and a reduced flux normalization, corresponding to an 
integral flux in the energy range (0.8 < E < 48.6 TeV) of 
(1.9+0.7)xl0- 13 cnr 2 s _1 . 



1 Obtained with t he data analysis pipeline CIAO v4.2 sherpa v2 
( Free man et all20 lj) : http ://cxc . cf a. harvard. edu/sherp a/. 
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Fig. 2. VHE y-ray spectrum of HESS J 1747-248 with 1 cr error 
bars, fitted with a power-law model. The fit results are discussed 
in the text. 



3. Discussion 

In the vicinity of the discovered VHE y-ray source, the archival 
multi-wavelength (MWL) data and astronomical catalogues do 
not show any typical candidate VHE y- ray emitter (e.g. an S NR, 
a young energetic pulsar or PWN, see IClapson et al.ll2.Ql lh . Its 
properties are unexpected for a GC, notably its extension, indi- 
cation of offset (~ 2.0o" level) from Terzan 5 and misalignment 
with the centre of the GC. Thus the source is challenging to in- 
terpret. 

Several approaches to constrain possible origins of the VHE 
radiation are considered here. First, chance coincidence prob- 
abilities, notably with young unrelated PWNe, are derived in 
Sect. 13.11 Next, the possible emission mechanisms for GCs are 
addressed. VHE leptons are covered in Sect. 13.21 in relation to 
msPSRs. A comparison with other GCs, identified as potential 
VHE y-ray sources, is performed in Sect. 13.31 Finally, possibili- 
ties for a hadronic emission scenario are discussed in Sect. 13.41 



3. 1 . Chance coincidence 

The proximity on the sky of the GC and the VHE y-ray source 
could simply be the chance coincidence of two physically un- 
related objects. In particular, the source parameters (extension, 
photon spectrum) are compatible with those of identified VHE 
PWNe (see e.g. iMattana et"ai1 l2Q09h . The probabilities of a 
chance coincidence based on the distribution of known objects 
are considered in two ways. 

The probability of source confusion in the VHE domain can 
be derived from t he distribution of sou rce counts as a function of 
Galactic latitude. fChaves et~al 1 (120091) find that this distribution, 
in the longitude range -85° < / < 60°, can be described by a 
Gaussian profile (containing 48 sources) centred at b — -0.26° 
with a width <x = 0.40° and additionally four outliers at latitude 
b < -2°, below the Galactic plane. There are no other detected 
VHE y-ray sources in the latitude band 1.5° < b < 2°. At a lati- 
tude of 1.7°, Terzan 5 is almost 5 <x away from the centre of the 
Gaussian, in a region where the total number of expected sources 
is neglig ible (< 10~ 3 ); i.e. the source presented here would be an 
extreme outlier in the Gaussian distribution. Using one source 



as the total expectation for this latitude band over the longitude 
range observed by the H.E.S.S. Galactic Plane Survey, the prob- 
ability of finding a VHE y-ray source within a conservative dis- 
tance of 0.1°from Terzan 5 is estimated to be ~ 10~ 4 . A caveat 
to this argument is that the sky coverage by H.E.S.S. is highly 
non-uniform, with longer exposure (favouring source detection) 
in regions typically located very close to the Galactic plane. 

A connection has been suggested between VHE y-ray 
sources and the PWN of pulsars with distance-scaled spin - 
down luminosity above 10 34 erg s _1 kpc 2 dCarrigan et al.l2008h . 
They could be observable during certain stages of their evo- 
lution with character istics similar to the detected VHE source 
(IMattana et al.l 120091) . In the band defined above, there is only 
one such pulsar (PSR J1 124-5916) in the ATNF catalogue 2 
dManchester et al]l2005l) . It has been proposed that the number 
of radio-loud and rad io-quiet pulsars might be comparable (e.g. 
iGonfhier et al.ll2004l) . therefore, one may expect a single unde- 
tected radio-quiet powerful pulsar. The resulting probability of 
chance coincidence with the VHE y-ray source within the rect- 
angular region shown in Fig. [T] is again ~ 10~ 4 . This approach 
is limited by the strength of the association (not every power- 
ful pulsar is associated with a VHE y-ray sour ce) and by the 
biase s and incompleteness of the catalogues (e.g. lGonthier et all 
12004 ), illustrated by the discovery of radio-quiet GeV y-ray pul- 
sars ( Abdo et aill2010bl but none so far in the neighbourhood of 
Terzan 5). 

In summary, it is unlikely that the association of the VHE 
y-ray source and the GC is by chance, but due to simplifications 
in the probability estimation, such a possibility cannot firmly be 
excluded. 



3.2. Leptonic VHE y-ray production 

Owing to the highly energetic electrons produced by their 
msPSR populations, GCs could be y-ray emitters, as predicted 
by three different models. In the GeV range, VJ09 explore the 
y-ray emission by the superposition o f pulsed rad i ation from 
individual pulsars, whereas BS07 and ICheng et alJ (1201 Oh ex- 
plore IC up-scattering by electrons either in the PWNe them- 
selves or elect rons that have bee n re-accelerated in colliding 
PWNe shocks. lAbdo et al.l d2009h find that the model of VJ09 
is in good agreement with the GeV properties of the Fermi- 
LAT source coincident with the GC 47 Tucanae, while the 
model of BS07 cannot rep roduce its spectral shape. More re- 
cently, ICheng et al.l (1201 Ol) have modelled the GeV spectra of 
47 Tucanae and Terzan 5 with IC emission. 

In the VHE range, these models rely on IC up-scattering 
of photons, from the cluster as well as the Galactic stellar 
and cosmic microwave backgrounds, either collectively by elec- 
trons leaving the msPSRs (see VJ09) or again by electrons re- 
accelerated in colliding shocks between PWNe (see BS07). The 
models of BS07 and VJ09 predict a flux in the VHE range of 
~ 1% of the Crab flux for reasonable input parameters, sim- 
ilar to the source fl ux determined in this work. For Terzan 5, 
ICheng et alJ (1201 Oh do not predict significant y-ray emission at 
energies beyond 100 GeV. Adopting the model of VJ09 with 
updated physical parameters for the GC dLanzoni et al.l [2~Q~Toh 
and a distance of 5.9 kpc, a population of 220+50 msPSRs for 
a magnetic field of 10 yuG (167+39 msPSRs for 20 ,uG) is re- 
quired to explain the observed luminosity of the VHE y-ray 
source above 1 TeV. This agrees w ell with the number of msP- 
SRs derived in I Abdo et al.1 ( 2010al) in the GeV range. The num- 



http://www.atnf.csiro.au/research/pulsar/psrcat/ (version 1 .40) 
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ber of required msPSRs is quite insensitive to the magnetic field 
in the range 10-20 fiG, although it increases rapidly outside 
this range (e.g. by ~ 50% at 5 fiG, see VJ09). This model ne- 
glects y-ray emission from outside r^, due to the uncertainties 
in the diffusion coefficient (and derived residence time) of the 
VHE electrons in this region, likely to be unusual if th e GC i s 
a dwarf galaxy remnant as proposed by iFerraro et al ] (I2009h . 
The additional interactions included in BS07 result in an ex- 
tended IC source of 2' - 3'radius, in rough agreement with 
the observed source size. When scaling their predictions for 
Terzan 5 ( their Fig. 7, solid curve , top panel) in the same way 
as done in lAharonian et all d2009l) . for 180 msPSRs of individ- 
ual spin-down power 10 34 ergs" 1 , the predicted integral flux is 
F(E>440GeV)^lxlO" 12 cm^s" 1 , at the level of the H.E.S.S. 
measurement. 

The IC emission shoul d be accompanied by synchrotron 
emission in the X-ray band (IVenter et al.ll2008h . Non-thermal X- 
ray emission e xtending beyond r^ has indeed been discovered 
from Terzan 5 feger et al.ll2010t) . The nature of this X-ray emis- 
sion, apparently centred on the GC, is still unknown and may 
be unrelated to HESS J1747-248. Based on the energetics a syn- 
chrotron origin seems to be preferred dEger et al.l 20101) . This 
interpretation is challenged by the spatial extension and the indi- 
cation of an offset of the VHE y-ray source from the GC centre. 
A simple model where the X-ray and VHE y-ray emission orig- 
inates in the same population of electrons cannot produce IC y- 
ray emission displaced from the peak of the stellar photon field 
at the GC centre as observed. 

To conclude, the observed flux - but not the morphology - of 
the VHE y-ray source is reasonably reproduced by the msPSRs 
scenarios. To explain the morphology, more sophisticated mod- 
els should be tested including particle trapping and non-uniform 
diffusion across the GC. Since the optical to near-infrared stellar 
photon field should be up-scattered by the VHE electrons, Klein- 
Nishina suppression of the IC process should become significant 
at multiple TeV energies. This would cause a steepening in the 
VHE y-ray spectrum, which should therefore not follow a pure 
power law. Statistics are currently too low to test this possibility. 



3.3. Comparison to other globular clusters 

Under the assumption that the observed VHE y-ray source is re- 
lated to Terzan 5, results from other GCs in the same energy band 
would be valuable. From the Fermi-LAT results (lAbdo et al.l 
l2010al) . NGC 6388 seems to host a number of msPSRs com- 
parable to Terzan 5. Its larger radiation field, even at its greater 
distance, makes it an intere sting candidate for VHE y-ray studies 
(lAbramowski et alT1l201 ll) . However, predictions of VHE y-ray 
fluxes for GCs strongly depend on uncertain parameters like the 
diffusion coefficient and the cluster magnetic field, as well as on 
the distance from the GC centre. 

Among the other GCs observed in the VHE domain, 
47 Tucanae is a very interesting object to compare with Terzan 5. 
It is located at 4 kpc (closer than Terzan 5) and away from 
the Galactic plane, thus simplifying observation at a ll wave- 
length s. It contains the second largest population (23, iRansoml 
120081) of individually detected msPSRs among GCs (after 
Terzan 5 with 33), in a greement w ith the Fermi-LAT esti- 
mate of 33+15 msPSRs (lAbdo et al. 



2010al). H.E.S.S. deri ved 



e.g. the distance to the GCs and the number of msPSRs, and the 
actual source morphology, this comparison only illustrates that 
additional observations of GCs will be required to settle the issue 
of VHE y-ray emission from their core. 



3.4. Hadronic VHE y-ray production 

The remnants of a type la supernova resulting from the merger of 
two white dwarfs could produce hadronic CRs (see e.g . the VHE 
y-ray detection of SN 1006 by lAcero etal] 1201 Obi). GCs are 



expected to boost the rate of stellar collisions ( Shara & Hurlevl 
l2002HGrindlav et al.ll2006l) . 

The energy in hadronic CRs required to explain a VHE y- 
ray source of luminosity L y with a hadronic scenario is E pp = 
LyTpp7]~ l (l + S), where r pp is the cooling time due to inelas- 
tic p-p collisions (linearly dependent on the density n of the 
ambient medium) and rj ss 1/3 is the fraction of CR energy 
that is converted into ji° mesons. CRs below the energy range 
probed by the observed VHE photons are accounted for by 
S . Here, the H.E.S.S. observation threshold of 440 GeV trans- 
lates into a minimum CR energy of 5 TeV. With n o indication 
of mo lecular material at the location of T erzan 5 (iDame et al. 
1 2001 ). n x 0.1 cm" 3 is assumed (following iDickev & Lockman 



a flux upper limit on 47 Tucanae dAharonian et al. I I2009T) of 
F(E>800GeV)<6.7xlO" 13 cm" 2 s" 1 , at 99% confidence level, for 
a point-like source. The flux from a point-like source analysis 
centred on Terzan 5 in the same energy range (see Sect. [2]) is be- 
low this limit. Given the uncertainties on the correcting factors, 



119921) . For a t ypical SNR CR spe ctrum below 5 TeV of spec- 
tral index 2.0. lAtovan et al.l d2006l) derive S « 5, giving E pp w 
10 51 (n/0.1cm -3 ) -1 (fif/5.9kpc) 2 erg, which is somewhat high for 
a supernova. 

An alternative scenario may be provided by the remnant of a 
short gamma-ray burst (GRB) in duced by the mer ger of two neu- 
tron st ars, see e.g. lNakarl(l2007l) for a review and Grindlav et al.l 
d2006h in the context of GCs. In relativistic shocks, the transfer 
of the initial kinetic energy to CRs is expected to be very efficient 
(lAtovan et al.l l2006). so the kinetic energy of a short GRB could 
roughly provide the energy in CRs needed for a hadronic inter- 
pretation of the H.E.S.S. source. Such a GRB remnant would be 
spatially extended, owing to the diffusion of CRs away from the 
explosion site, in relation to its age t « r 2 /2D, where r is the 
radius of the remnant and D the diffusion coefficient. 

Assuming CRs of energy 5 TeV, the extension of the present 
source would give t « 10 3 yr fo r D = 10 28 cm 2 s"' (in the 
Galactic disk, lAtovan et aT1l2006l) . The rate of short GRBs in 
the Galaxy should be about one event every 10 4 (/ b "'/100)" 1 yr, 
where /b is the GRB beaming fa ctor, uncerta in but expected to 
fall in the range 1«: <100 dNakarll2007l) . This estimate is 
based on a rate of short GRBs of 10 Gpc" 3 yr" 1 (lNakarll2007l) 
and a density of Milky Way-ty pe galaxies in the local Universe 
of 10" 2 galaxies per Mpc" 3 dCole et al.l 1200 lh . For D below 
10 27 cm 2 s" 1 in the TeV range in the vic inity of the GC, as sug- 
gested by ICrocker & Aharonianl (1201 Ol) . the age obtained from 
the source extension is roughly compatible with the rate of short 
GRBs in the Galaxy. Nevertheless, the value of D in this region 
may differ significantly from the available estimates. 

A hadronic scenario could accommodate the observed pure 
power-law spectrum in the VHE y-ray domain, as opposed to a 
leptonic scenario. Detection of thermal X-rays would strengthen 
the case for a GRB remnant, where shocks driven by sub- 
relativistic ejecta e xpelled during the merger would heat the in- 
terstellar medium dDomainko & Ruffertll2005l 120081) . However, 
no observational evidence was found in support of either the su- 
pernova or GRB remnant scenario. 
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4. Outlook 

The nature of HESS J1747-248 is uncertain, since no counter- 
part or model can fully explain the observed morphology. An 
association with Terzan 5 is tantalizing, but the available data do 
not firmly prove this scenario. 

Several tests could be done to clarify the nature of the VHE 
source. In the X-ray range, a large FoV and an extended en- 
ergy range would be desirable. XMM-Newton with its large FoV 
is suitable for searching for counterparts over the whole extent 
of HESS J1747-248. As a com plement, the upcoming NuSTAR 
mission (lHarrison et al. 1 1201 Ol) will provide significantly im- 
proved sensitivity in the ~ 10-70 keV range, which is less 
affected by neutral hydrogen absorption. This may help in con- 
straining the nature of the diffuse X-ray emission found by 
Chandra. Pulsar searches, with Fermz-LAT and radio telescopes, 
may reveal a powerful pulsar in the vicinity of HESS J1747-248. 
Finally H.E.S.S phase II, which will extend observations in the 
VHE y-ray range to lower energies, will further characterize the 
source. 

In parallel, more sophisticated source models might narrow 
down the list of applicable emission scenarios. 
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